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A model of inflation is presented where the inflaton field is a complex scalar field coupled to 
\^ \ a C/(l) gauge field. Due to the axial symmetry of the potential, the inflation is driven by the 

radial direction while the angular field is gauged by U{1). Due to the coupling of the inflaton 
to the gauge field, a time dependent mass term for the gauge field is generated dynamically and 
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, conformal invariance is broken. We study whether a significant amount of primordial magnetic 
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fields can be generated during infiation by allowing a time-dependent U{1) gauge kinetic coupling. 
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I. INTRODUCTION 



Astronomical observations indicate the existence of magnetic fields in galaxies and also 
in cosmological scales with coherence lengths as big as 1 Mpc and with the magnitude of 
order 10^^ Gauss, for more details see l|-l3|]. The origin of the cosmic magnetic fields are 
not well-understood. One possible explanation is that they may have primordial origin 
which are later amplified by galactic dynamo mechanism into the current observed value. 
The estimation for the magnitude of the primordial seed is not certain. Conventionally, it 
is assumed that a seed in the range 10~^^ — 10~^^ Gauss at the time of matter-radiation 
decoupling may be required for dynamo mechanism to produce the current observational 
value. This lower bound may be relaxed to 10^^^ Gauss in a flat dark energy dominated 
Universe j4|. ^ 

Starting with the work of Turner and Widrow [Sj], during past two decades there were 
some attempts to obtain primordial magnetic seed through inflation j6-43|. In particular, 
inflation can naturally stretch small micro-physical scales inside a casual patch into the 
cosmological sizes. This can provide an elegant mechanism to explain the coherence of mag- 
netic fields on distances much bigger than astrophysical scales. The key obstacle, however, 
in producing magnetic field during inflation (in general in a Friedmann- Robert son- Walker 
(FRW) background) is the conformal invariance. To produce magnetic field one has to break 
the conformal invariance of the classical electrodynamics. There are several mechanisms to 
break conformal invariance. These include: coupling electrodynamics non-minimally to grav- 
ity which produce a time-de penden t mass for photon js], [igI, introducing a time-dependent 
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kinetic coupling 
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43| , and couphng the photon to a charged scalar field 



16|. 



In this paper we combine the last two mentioned methods together in order to break the 
conformal invariance and obtain an appreciable amount of seed magnetic field. We consider 
the model where the inflaton field is a complex scalar field coupled to a U{1) gauge field. 
The conformal invariance is broken dynamically when the scalar field is displaced from its 
minimum during inflation. Furthermore, the time-dependence of the gauge kinetic coupling 
as in 



22 



43| can help to produce a considerable amount of primordial magnetic fields seed 
required for galactic dynamo. As indicated recently in j43|, the key constraints in models 
with time- varying gauge coupling are (a): the requirement that the time-dependent gauge 
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coupling remains small for all time during inflation, so the perturbative gauge theory is 
under control and (b): one must check that the strong back reaction from gauge fleld does 
not destroy the inflationary background. We shall see that these constraints in addition with 
the requirement to have sufficient inflation basically control the amplitude of the magnetic 
flelds seed produced. 

This paper is organized as follows. In section [TTl we present our set up. In section UTTl the 
magnetic flelds produced in our model are studied followed by discussions and conclusion in 
section IIVI While this paper was in its flnal stage, the work [44,] appeared which has some 
overlaps with our studies. 



II. SET UP 



Here we present our set up. It contains a complex scalar fleld, 0, coupled to a U{1) gauge 
fleld, Af^, with a time dependent gauge kinetic coupling, I(t). The action is 



s 



P 1 P(f) 



(1) 



where Mp^ = SnG and G being the Newton constant. The covariant derivative is given by 



D^4> = d^cf) + ie(j)A 



^l 1 



(2) 



where e is the dimensionless gauge coupling of to 0. As usual, the gauge fleld strength 
is given by 



duAy — dyAa ■ 



(3) 



Motivated by 22|, |43| we have entertained the possibility that the gauge kinetic coupling is 
time dependent. This may arise for example in string theory setup where the gau ge kinetic 
action is coupled to the dilaton which is running with time. As suggested in 43| one may 
also consider that the time dependence of I{t) originates from the time dependence of e{t). 
This may be obtained by replacing A^ — )■ A^/e{t). But this has the problem that the kinetic 
action F^^F^^^ becomes modifled under A^ — A^/e{t) and the gauge invariant is lost. At 
the phenomenological level, we shall take e to be constant and assume the time dependence 
of I{t) originates from the coupling of the gauge kinetic action to other flelds such as dilaton. 
However, the dynamics of the dilaton should be such that it does not destroy the slow-roll 



properties of the inflationary pot ential. This may be a non-trivial task, but we shall take 



the phenomenological view of 
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43| and proceed with action represented in Eq. ([T]). 



We work with potentials which have axial symmetry where V is only a function of < 
p. It is more instructive to decompose the inflaton field into the radial and angular parts 



(4) 



so V = V{p). As usual, the action, Eq. ([T]), is invariant under local gauge transformation 



^^^^ 



e e + e{x). 
With this decomposition, Eq. ([T]) is transformed into: 
S = I d'^Xy/—g 



The corresponding Klein-Gordon equations of motion are: 





= [V^d^p] 



J.J'' 
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accompanied by with the Maxwell's equation 

d,{^gI{tfF''^)=er 
where the current, J^, is defined by 



r = p'^^d'^e + eA' 



(5) 

(6) 

(7) 
(8) 

(9) 
(10) 



The conservation of from Eq. ([7]) is a manifestation of the axial symmetry imposed on 
V . Interestingly, Eq. ([7]) is not independent from Maxwell's equation, where F^'^ being 
anti-symmetic leads to d^di,F^^^ = 9^J^ = 0. 

Finally, the stress energy momentum tensor, T^/s, for the Einstein equation. Gap = 



SnGTnB, is: 



-I{tf 



ga^F^uF'^" + litfFa^^F^ + dapdpp + 



JaJl 



p'^\g\ 



9aP 



V 



'IV. 



At the background level, we start with the isotropic and homogeneous FRW space-time 
with the metric 



-df + a{tfdx 



(12) 



Considering the u = component of the Maxwell equation (Q implies that 

jO = pV(^ + eA°) = 0. (13) 



This indicates that the total electric charge is zero [45|. This is somewhat similar to usual 
angular momentum conservation where in the absence of gauge field, an axial symmetric 
potential leads to angular momentum conservation dt{a^p'^6) = 0. 

The gauge field is invariant under gauge transformation Eq. ([5]). We need to fix 
the gauge to study the physical independent degrees of freedom. We choose the Coloumb- 
radiation gauge where Aq = d^Ai = 0. With Aq = 0, from J° = 0, one obtains that 6 = 
at the level of background. The z/ = component of the perturbed Maxwell equation (Eq. 
([9])) with the above gauge implies that 69 = 0. This is also consistent with the perturbed 
angular field equation (Eq. ([7])). 

A non-zero spatial component of vector field at the background level may produce large 



anisotropies. In [4J] the mechanism of magnetic seed production in the presence of an 
anisotropic gauge field was studied to some extent. It was shown that an anisotropic gauge 
field background will modify the inflation dynamics. It is found that the final produced 
magnetic field becomes highly suppressed as compared to the isotropic background. Also it 
is argued that a large anisotropic background may produce large anisotropies on CMB which 
may be detectable and need further studies. Here, in order to simplify the analysis and. 



in the light of 4^, to have an optimum primordial magnetic field production, we assume 
that A = at the background level. We need to check if this is a consistent solution of the 
background inflationary dynamics. To this end, we note that in our isotropic background, 
Fij = J* = and one can easily check that A = A = is a consistent solution of Eqs. 
([T])-® and the Einstein equations. In other words, the p field and the gauge field equations 
of motion decouple from each other. One also has to check that this ansatz is stable against 
perturbations. To see this, we note that the mass squared term created for the gauge field 
via spontaneous symmetry breaking (the Higgs mechanism) is positive, given by 
which indicates there is no tachyon in the mass spectra of the gauge field. 

With this ansatz, at the background level, the independent equations of motion are 

= p + 3Hp + Vp 
3H'Mj, = \(? + V{p) . (14) 
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The interesting result is that the inflation is completely driven by the radial fleld p and the 
angular fleld 6 is gauged by U (1). As in conventional single fleld models of inflation, one can 
choose V{p) to be flat enough to support an extended period of inflation. To be specific, 
below we consider the chaotic inflationary models. For potential V = the number 

of e-folding, A^e, is given by 4 Mj, Ne — pl-, where pi is the initial value of the scalar fleld. To 
solve the flatness and horizon problem, we may take Ne = 60, corresponding to pi ~ lOMp. 
Furthermore, to flt the WMAP normalization for the density perturbation, P-ji ~ 2.4 x 10~^, 
one requires that m ~ 6 x lO^^Mp. For the quartic inflationary potential V = A0^/4, to flt 
the the data one requires that A ~ 10^^'^ and pi ~ lOMp. 



III. MAGNETIC FIELD PRODUCTION DURING INFLATION 



Here we study magnetic fleld production during inflation. As is well-known, magnetic 
fleld production in FRW backgrounds is heavily suppressed due to the conformal invariance. 
To produce a signiflcant amount of magnetic fleld, conformal invariance of the Maxwell 
equation should be broken. As we see from Eqs. ([6]) and (Q, the conformal invariance is 
broken via spontaneous breaking of U{1) symmetry where a mass term for the photon is 
generated dynamically through the coupling of p fleld to A^. The possibility of generating 
magnetic fleld during inflation due to t/Yl) symmetry breaking by a complex fleld is studied 
extensively in literature , |l8| . In these models the complex scalar fleld is usually 



different than the inflaton fleld and it is assumed that the mechanism of magnetic fleld 
production happens on top of the inflationary background without affecting the dynamics 
of inflation. However, in our model we assume that the inflaton fleld is the same as the 
complex scalar fleld, i.e. the fleld p. This has the advantage that the parameters controlling 
the magnitude of magnetic flelds production are directly related to the parameters that 
control the inflationary predictions, such as N^. and density perturbations. 

The vector perturbations are generated only quantum mechanically during inflation. To 
study their evolution and production, we go to Fourier space where the equation for the 
spatial components of the vector fleld is 

A'!, + 2jA[, + {k' + Y2 A,k = ^, (15) 
where k is the comoving wave number. Hereafter, the prime indicates the derivative with 
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respect to the conformal time dr = dt/a{t). This equation represents a damped harmonic 
oscillator with a time- dependent mass. One can get rid of the first derivative term in Eq. 
([IS]) using Vk = /(r)i'fc, so 

v:+(k'-^-^ + ^^a'p'y, = 0, (16) 



For convenience we omits the indices i. As in 
gauge kinetic coupling 



43|] we consider the following ansatz for the 



Hr) = If(-] , (17) 



where If and a/ are the values of /(r) and the scale factor at the end of inflation, respectively. 
We expect that < 1. Noting that the gauge kinetic coupling is inversely related to /, a 
negative value of p corresponds to the case when gauge coupling is very large at the beginning 
of inflation and perturbative analysis is not reliable in the beginning of inflation. On the 
other hand, for a positive value of p the gauge coupling is very small at the beginning of 
inflation and with the assumption // ^ 1 the perturbative gauge theory is under control for 
all time. For this purpose, we take p > in our analysis below. Plugging the ansatz Eq. 
(II7D into Eq. (JISD yields 

where 



<+{k'-'-'^ + -^,]v, = 0, (18) 



/3 . (19) 

To obtain Eq. ( ITSl) the relation a ~ —1/Ht is used during the slow-roll inflation, where 
H is the Hubble constant during inflation. Also we assume that p is changing very slowly 
during inflation so we replace by its average value (p^). 



The amplitude of the magnetic field , 6b, at the end of inflation is given by 43 1 



where \vk{Tf)\ is the magnitude of Vk at the end of inflation. 

Eq. ( ITSl) can not be analytically solved for arbitrary value of p. Before considering general 
positive value of p, first we consider the particular case of p = corresponding to a constant 
gauge kinetic coupling, / = 1. 



A. Constant Gauge Kinetic Coupling, p = 



For p = the solution of Eq. (fT8|) are given in terms of Hankel functions 

Vk = Ak 



2 

where z/q = 1/4 — /3. To match the initial vacuum state 



(21) 



"' = 71f '''' 



one requires that bi = 1 and 62 = in Eq. ( 1211) . On the other hand, for long wavelength 
modes, k\T\ — )• 0, and by using the asymptotic relation 



HS^k\r\)^-^^{uo)(^] , (23) 



at the end of inflation one obtains 



r(z/o) 2^^0-3/2 fafH 



1/0-5/2 



where T{x) is the Gamma function. Noting that the real value of uq can not exceed 1/4 the 
amplitude of magnetic field in Eq. is hugely scale suppressed. 

To calculate the amplitude of magnetic fields after inflation we should take into account 
the expansion of the Universe which dilutes the magnetic field energy density, -Bj-B*, like ra- 
diation and 6b decreases like a~^. Assuming that the preheating and reheating happens in- 



stantly, followed by a radiation-dominated Universe, one obtains that adec/o-f — VHMp/Tdec 
where Tdec and adec are the values of the temperature and scale factor at the time of de- 
coupling, tdec, respectively. For H ~ lO^^Mp one obtains a/ ~ 10^^^ and adec/ (if — 10^^ 
where we have set the magnitude of the scale factor today equal to unity, gq = 1. Putting 
all together, with ^ lO^^^Mj, ~ 10"^^ Gauss, the amplitude of the magnetic fields at 
the time of decoupling is 6B{tdec) ^ 10'^ X 1011(2^0-5) < ^0^50 Gauss on comoving scales 
= IMpc ~ 10^^GeV~^. This seed magnetic field is too small to be amplified by galactic 
dynamo mechanism into the current observed values. 

One may compare the case p = here to the models studied in [sl, Q, 43 1 where the 



gauge field is coupled to the Ricci scalar via the mass term l3RAf^A^ /2 and the gauge 
invariance and the conformal invariance are broken explicitly. However, in our model the 
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gauge invariance and the conformal invariance are broken only dynamically (spontaneously) 
when the infiaton field is dislocated from the minimum during infiation. In their cases the 
conformal coupling parameter /3 can have a negative value so a negative mass squared can 
be generated for the gauge field. This results in an amplification of the magnetic field, and 
one has to check that its strong back reactions do not destroy the background infiationary 
setup [43|. However, in our case f3 > and a positive mass squared term is created for the 
gauge field and the magnetic field is more suppressed compared to the case of a massless 
gauge field where (3 = and z/ = 1/2. 

B. Magnetic Fields from Time Varying Gauge Kinetic Coupling, p > 

To keep the gauge kinetic coupling under perturbative control, we take p > 0. As it is 
obvious, Eq. ( IT8|) can not be solved analytically, so we rely on the numerical solution and 
asymptotic behaviors. 

During early stage of inflation the first two terms inside the bracket of Eq. ( ITSll dominate. 
The last term dominates only at the final stage of infiation when r — 0. For the last term 
to dominate over the second term one requires that 

n^P\r,r^='-j0>>l. (25) 

To have an estimate of Q, consider the infiationary potential V = m?p^ /2. Using the relation 
GH'^Mp = m^p^, one obtains 

To fit the COBE normalization of density perturbation, m ~ lO^^Mp so, for e ~ // ~ 1, 
VL can be as large as 10^^. Increasing(decreasing) the ratio e/I{Tf) the magnitude of fl 
increases(decreases). On the other hand, for A0^/4 infiationary model, one obtains 

For A ^ 10 and p ~ lOMp, Q can be as large as 10^^ like in the previous example. 

In the early stage of infiation when the term containing /3 is negligible in Eq. f[T5|) . the 
asymptotic solution is 



where u = p — 1/2. At the final stage of infiation, the term containing /3 dominates and the 
asymptotic solution of Eq. f lTSj) is given by: 

^VW\ L ^(1) f v^l 



final / 



r 



Cl 



P 



(29) 



where /i = l/2p and Ci and C2 are the Bogoliubov coefficients. One needs to match solutions 



vl^"^^^ and ff"^' at the transition time, t^, where 



(3\rc 



-2p 



1)1 • 



(30) 



Since Tc is expected to be close to the end of infiation, then kr^ — )■ and given above is 
independent of k. 

Demanding that both Vk and v'f. to be continuous at fixes the coefficients Ci and C2 

±71 



Cl,2 



4p 



^i(l+2i/)7r/4 



(2,1) 



(31) 



where the arguments of if^ are a/^| rd ^/p while that of Hi, are \kTc\. 

It is instructive to introduce the dimensionless variable x = k\T\, and rewrite Eq. ( I30ll as 



l/2p 



Xf. 



(32) 



Jp(p- 1) 
where Xf = k\Tf\. 

Figure [U shows the value of Xc/xf as a function of p. For a given value of p, by increasing 
Q the value of Xc goes to the larger values as expected. We note that with k~^ = IMpc and 
H ^ lO^^Mp one has Xf ^ 10~^^. This indicates that Xc — ?■ as verified from Eq. fl30l) and 
also from Figure 1. Physically, this means that the last term in Eq. f[T8|) dominates at the 
late stages of infiation. 

The full Numerical solution of Eq. ( IT8|) for typical value of p = 2.3 and Q = 10^^ along 
with the asymptotic behaviors at the initial stage, Eq. (125]) . and the asymptotic solution 
at the final stage, Eq. ( !29|) . are plotted in Figure [21 One can check that the approximate 
solution Eq. (1291) is in a good agreement with the full numerical solution. The transition 
point, Xc, where the two solutions Eq. (!29|) and (!28|) are matched smoothly together can be 
seen in Figure 2. 

As explained above, we expect that kTc — 0, which can be used to simplify the value of 
Cl and C2. Using the asymptotic form Eq. ( I23l) . one obtains that 



Cl 



C2 



\k 



(33) 
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FIG. 1: Here we plot the behavior of Xc = /c|tc| versus p given by Eq. (f32l) for various values of 
17. We took Xf = 10-22. 



and 



.final \ 



(34) 



where to get the final answer use was made of the large argument limit of the Hankel 
functions H^{x) ~ x~^/'^ for x ^ 1. One the other hand, using Eq. ( l30l) one can express 
Tc in favor of Tf via Tc/tj ~ 0}/"^"^. Putting all together, and noting that v = p — 1/2 and 
neglecting factors of order unity one has 



Tj\VL 4 + 4p|/;;7- 



(35) 



which from Eq. fl20|) results in 



k 



p-3 



(36) 



As expected, up to term containing the factor of Q, this is the same as in 43] . Since 
^ 1, one rnay expect that our value of ^^(rdec) would be more suppressed compared to 



the result of 



43|. On the other hand, due to the presence of Q, the allowed range of p 
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FIG. 2: The numerical solution of the real part of Eq. (llSp as a function of x = A;|t|. To make the 
presentation more transparent, the vertical axis is scaled down by a factor of 10~^^. The solid (red) 
curve indicates the full numerical solution whereas the other two curves indicate the early and the 
late time solutions, given by Eqs. ([28]) and ([29]). We took p = 2.3, n = 10^^ and x/ = 10-^2 . The 
transition point, Xc, calculated from Eq. (pO|) . is indicated by the vertical line. We see that at Xc 
the solutions (|28p and (j29p match smoothly and are in a good agreement with the full numerical 
solution. 



increases as compared to that in 43| . One has to check whether the increase in range of p 
can compensate the suppression from large Q effect. In the examples below we will show 
that unfortunately Eq. f l36|) gives smaller value for the magnetic fields than those obtained 



m 



43 1 in the absence of fi. 

n 

As in [43i|, there is an upper bound on p, determined by the requirement that the energy 
density from the vector field, eEM, does not exceed the background inflationary energy 
density. This is equivalent to e em < Mp H"^. There are two sources for eem coming from Eq. 
ffTTj) . The first one, which was also considered in |43|, is ejf^/ = -\l{tfF^pF''^ + I{tfFQ^F^ 



and the second one is the contribution from the gauge coupling e in the form of e^}^ = 
JaJ"" / p'^\g\ = Below we show that e^^\j- dominates over by a factor of Q. 
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Following 43|, one obtains 



-EM 



(rj) = -eV(0|A,A^|0) 



fH 



dke\vk{rf)\' 



(37) 



where at is the initial size of the scale factor at the start of inflation. Using Eq. ( IMI) . one 
obtains 



El 

4^ 



X < 



1 

2-p 

In 

1 "/ 
p-2 \ aj 



2{p-2) 



p < 2 
p = 2 

p > 2 



(38) 



l/ajH is used during slow-roll inflation. 

arH 



where to get the final result, the relation Tf ^ 

On the other hand, is given by 43| ~ aj^ XTiT '^^ I'^U''"/)^- show that 

e~^li ~ l^^^e^jj^^ and with 1, the requirement e~^j^j < MpH"^ is the stronger constraint 
that has to be satisfied. 

From Eq. ( l36l) . to have a large enough magnetic field, it is optimum to have p > 2. 
However, the magnitude of p is controlled by the value of Q and af/ai. For example, 
with H ~ 10'^ Mp, Ne = ln(a//aj) = 60 and Q = 10^^ which is naturally obtained in 
chaotic inflationary models, the constraint e^^\j < MpH"^ results in p < 2.3. This results in 
^Bitdec) ~ 10~^^ Gauss on ~ Mpc scales. This value of magnetic seed may be marginally 
acceptable in the light of {4]. Reducing results in an increase in ^^(tdec)- For example, 
choosing Q = 10^, one obtains 6B{tdec) ~ 10^^^ Gauss. With // ~ 1, this value of fl 
corresponds to e ~ 10~^ which indicates a severe fine-tuning. On the other hand, in the 
absence of Q factor, one obtains p < 2.2 and 5_B(tdec) ~ 10^^^ Gauss. 



IV. CONCLUSION 

We considered a model where the inflaton field is a charged scalar field coupled to a 
U{1) gauge field with a time-varying gauge kinetic coupling. The gauge field is frozen in 
the classical background and is excited quantum mechanically. The amplitudes of magnetic 
fields produced down to the end of inflation and at the time of decoupling are calculated. 
Our model has features in common with models such as in jsf 
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14j . IigI . [isl where a charged 



scalar field is present in an inflationary background. However, since in our model the charged 
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scalar field itself is the inflaton field, the prediction of magnetic fields are directly controlled 
by the inflationary parameters such as e, A^^e and m in m^p^ or A in X(f)'^ models. As in 43| . 
the main constraint controlling the magnitude of 5^ comes from the requirement that the 
back reaction from the gauge field does not destroy the inflationary background. This in 
turn imposes an upper bound on the range of parameter p. Our model predicts values of 
^sitdec) somewhat smaller than what is obtained in |43|] in the absence of charged scalar 
field. However, our model predicts a higher value of spectral index parameter, p, for the 
primordial magnetic field. With natural parameter values, our model predicts that magnetic 
field at the order of ^^(tdec) ~ 10^^^ Gauss can be created on k"^ ~ Mpc scales. This value 
may be marginally acceptable in the case of very efficient dynamo mechanism. 

In this model important issues such as the amplification of magnetic fields during pre- 
heating and subsequent suppressions via electric conductance are not considered. Specially, 
since the gauge field obtains a time-dependent mass, as is evident in Eq. (fT5l) . the effects 
from parametric resonance and significant amplification of magnetic field during preheating 
stage may play some important roles 1^, 14, 16]. This may help to relax the bounds above 



by few orders of magnitude. 

In this work we have turned off the background gauge field classically. Although this 
is a consistent solution, but it is not the most natural solution and requires an additional 



fine-tuning in our model. In a future work 



461] we study the case in our model where the 



background gauge field is not turned off classically. For this, one has to search for the 
parameter space where the gauge-field energy density does not dominate over the inflaton 
field energy density such that the produced anisotropics are within the observational bounds 
from CMB. This is similar to the analysis of jl^] except that in their model there is no gauge 
coupling, e = 0, and /3 = 0. 
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